A new method for numerical differentiation based on direct and inverse problems of partial differential equations

In this paper, we mainly study a numerical differentiation problem which aims to approximate the second order derivative of a single variable function from its noise data. By transforming the problem into a combination of direct and inverse problems of partial differential equations (heat conduction equations), a new method that we call the PDEs-based numerical differentiation method is proposed. By means of the finite element method and the Tikhonov regularization, implementations of the proposed PDEs-based method are presented with a posterior strategy for choosing regularization parameters. Numerical results show that the PDEs-based numerical differentiation method is highly feasible and stable with respect to data noise.     

Numerical study on electrostatic potential distribution of large cylindrical oil tanks

This study carried out a numerical simulation of electrostatic potential distribution in the cylindrical oil tank, focusing on spatial distribution of electrostatic potential and the variation tendency of the potential in the filling operation. Electrostatic potential distribution of stationary oil was illustrated, and the relationship between the electrostatic potential distribution and some parameters of the tank and oil was analyzed as well. Dynamic electrostatic distribution was studied when charged oil is filling into the tank, and ESD prevention measures’ effects on electrostatic potential distribution were simulated. Finally, calculation results of two other types of oil tanks were presented in brief.     

力学问题加权残值法的双侧优化逼近

本文给出的力学问题加权残值法的数值解法,可直观地表示出解的最小上边界和最大下边界。     

采动诱发岩体移动破坏过程数值模拟研究

本文基于岩石介质的宏观非线性主要是由非均质性和各向异性造成的, 应用新的数值计算软件RFPA(2D), 对采动引起岩体失稳破坏的全过程进了数值模拟研究。     

宝泉抽水蓄能电站岩体抗剪参数的选取

本文讨论了岩体抗剪参数C、f值的选取方法, 推导了C、f值的数字特征公式。通过对宝泉抽水蓄能电站工程岩体抗剪参数的选取, 结果表明：由于岩体抗剪参数具有随机性和模糊性, 采用随机-模糊法选取的C、f值要符合实际；文中给出的数字特征公式可对C、f值的波动性及相关性进行检验。     

Compact finite difference-Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported. The project supported by the National Natural Science Foundation of China     

两斜交圆柱壳应力分析的边界元法

本文研究了两斜交圆柱壳的边界元法。提出了用板的基本解叠加级数形式的修正项构成圆柱壳的基本解,提高了计算精度,缩短了计算时间。对区域积分进行特殊处理,从而避免了内部网格的划分,大大减少了数据准备工作量和占机内存。并编制了FORTRAN计算程序,进行了数值计算。     

样条边界元法分析圆形基础板非线性非光滑接触

本文用样条边界元法分析轴对称圆形基础板的非线性非光滑接触。除事先给出外力与地基表面位移的对应关系外,地基的模型、板的厚薄(Reissner/Kirchhoff型)和外荷载沿径向的分布均可是任意的。文中将地基反力看作板的外载,对变量沿径向采用不等距祥条插值,并用无地基效应的板基本解。算例表明本文方法只要用少量自由度,便可得到令人满意的结果。     

Parametric identification of nonlinear systems using multiple trials

It is observed that the harmonic balance (HB) method of parametric identification of nonlinear system may not give right identification results for a single test data. A multiple-trial HB scheme is suggested to obtain improved results in the identification, compared with a single sample test. Several independent tests are conducted by subjecting the system to a range of harmonic excitations. The individual data sets are combined to obtain the matrix for inversion. This leads to the mean square error minimization of the entire set of periodic orbits. It is shown that the combination of independent test data gives correct results even in the case where the individual data sets give wrong results.